library(ggplot2)
1:16
x <- 1:16
y <- 1:16
df = data.frame(list(x, y))
df
ggplot(data = df, aes(x = factor(x), y = y, shape = factor(x))) + geom_point() +
    scale_shape_manual(values=factor(x))

